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Bilayer graphene nanoribbon with zigzag edge is investigated with the tight binding model. Two 
stacking structures, a and f3, are considered. The band splitting is seen in the a structure, while 
the splitting in the wave number direction is found in the ft structure. The local density of states 
in the ft structure tend to avoid sites where inter-layer hopping interactions are present. 



I. INTRODUCTION 

The graphite, multi-layer, and single-layer graphene 
materials have been studied intensively, since the electric 
field effect has been found in atomically thin graphene 
filmsii These materials can be regarded as bulk systems. 
On the other hand, nanographenes with controlled edge 
structures have been predicted to have localized states 
along the zigzag edges^ The presence of the edge states 
has been observed by experiments of scanning tunnel- 
ing spectroscopy^ Thus, the studies of the edge states 
are one of the interesting topic of the field. The recent 
atomic bottom-up fabrication of nanoribbons really pro- 
motes experimental and theoretical investigations. 5 

Previously, one of the present authors has studied 
the stacking effects of the nanographene by considering 
weak inter-layer hopping interactions in the tight binding 
model*^ The cluster calculations have been performed, 
and compared with experiments of the magnetic prop- 
erties. It has been found that open shell nature of the 
electronic orbital of the each layer is important. The 
existence of the spin in the each layer gives rise to the 
magnetism of the stacked nanographene. Therefore, the 
simple tight binding model is effective. 

In this paper, bilayer graphene nanoribbon with zigzag 
edge is investigated with the tight binding model. Two 
stacking structures, a and /J, are considered. The band 
splitting is seen in the a structure, while the splitting in 
the wave number direction is found in the /3 structure. 
The local density of states in the j3 structure tend to avoid 
sites where inter-layer hopping interactions are present. 

This paper is organized as follows. In the next section, 
the model is explained. The results are shown in Section 
3. The paper is closed with summary in Section 4. 



II. MODEL 

We consider the following tight binding model 

H = - t E (4,^+h.c.) 




FIG. 1. A-B stacked bilayer graphene nanoribbon with zigzag 
edges. The upper layer is shown by the solid lines, and the 
lower layer by dotted lines. In the a structure (a), the upper 
layer is shift by the bond length downward to the position of 
the lower layer. The region surrounded by the dashed line is 
the unit cell in the direction of the one dimensional direction. 
At the circles, two carbon atoms of the upper and lower layers 
overlap completely, and there is the weak hopping interaction 
ti here. In (b), the j3 structure is shown, where the lower layer 
is shift right-down direction so the stacking pattern is different 
from that of the a structure. The convention of the a and ft 
structures are the same used in the paper.— (c) The detailed 
view around the edge atom in the /3 structure is displayed. 

-h J2 (4,<a*+h.c), (i) 

where Ci jCT and c?i j(T are the annihilation operators of elec- 
trons at the lattice site i of the spin a on the upper 
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FIG. 2. Energy band structures for the single layer (a), a 
structure (b), and /3 structure (c). Details near the Brillouin 
zone edge are magnified for the single layer (d), a structure 
(e), and ft structure (f). In (e) and (f), the energy bands of 
the single layer are shown by the red lines for comparison. 
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FIG. 3. Number of states against the energy for the single 
layer (a), a structure (b), and ft structure (c). 

and lower layers, respectively. The quantity t is the hop- 
ping integral of tt electrons between neighboring lattice 
sites. Two stacking patterns, shown in Figs. 1 (a) and 
(b), are considered. They are named as a and ft struc- 
tures, respectively. This convention has been used in the 
literature^ There are N z = 6 zigzag lines in upper and 
lower layers. The detailed view near the edge atom in the 
ft structure is displayed in Fig. 1 (c). There is a weak 
hopping integral t\ along the dashed line. 



is ti = 0.1/;. In (b), the split of the energy bands is 
seen compared with Fig. (a). In contrast, split in the 
perpendicular direction is small in Fig. (c). In order to 
see split structures clearly, details near the Brillouin zone 
edge are magnified, for the single layer (d), a structure 
(e), and ft structure (f). In (e) and (f), the energy bands 
of the single layer are shown by the red lines for com- 
parison. The nearly flat band due to the edge state in 
— 7r < kd < — 27r/3 is present at the energy E ~ in (d), 
where d is the unit cell length of the one dimensional di- 
rection of Figs. 1 (a) and (b). The energy bands starts at 
E = l.Ot at kd — tt, typical to the graphene structure. In 
(e), the energy split is magnified again. However, in (f), 
the energy split is not seen, and split in the wave number 
direction is found. This property is confirmed by looking 
at the numerical data, also. 

It is interesting to look at how such the difference of 
the band split in the a and ft structures appear, in the 
quantization of conductance. Number of states in the en- 
ergy window for the positive wave number is calculated 
for the single layer, a, and ft structures. Fig. 3 displays 
the calculated results for the parameters N z = 20 and 
ti = O.lt again. Around the energy E = of the single 
layer (a), the energy band is singly degenerate, as shown 
in Fig. 2 (d). At the energy E — 0.4i, the energy win- 
dow reaches the first parabolic band, and the number of 
states jumps by two. The energy window reaches the 
second parabolic band at the energy E = 0.65i. In this 
way, a series of jumps by two is realized in (a). For the 
a structure (b), number of states near the energy E = 
is two due to the bilayer property. There are two steps 
of the jump of two around the energy E = 0.4i, owing to 
the energy split of the first parabola. Also, two steps are 
present closely around the energy E — Q.7t due to the 
second parabolic band. On the other hand, each jump of 
the number of states becomes four for the ft structure (c) 
due to the split in the wave number direction. The de- 
pendence of the energy reflects the band structures, and 
this will appear quantization of conductance experimen- 
tally. 



IV. SUMMARY 



III. RESULTS 

In order to see effects of inter-layer interactions, the 
tight binding model is solved numerically. The model is 
solved using the Bloch theory of the band calculation. 
Energy band structures of the N z = 20 systems are dis- 
played in Fig. 2, for the single layer (a), a structure (b), 
and ft structure (c). The interlayer interaction strength 



In summary, weak inter-layer interactions have been 
considered for the bilayer graphene nanoribbon with 
zigzag edge. The a and ft stacking structures have been 
considered. The band splitting is seen in the a structure, 
while the splitting in the wave number direction is found 
in the ft structure. The local density of states in the ft 
structure tend to avoid sites where inter-layer hopping 
interactions exist. 
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